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Abstract. - We investigate the single-point velocity probability density function (PDF) in three- 
dimensional fully developed homogeneous isotropic turbulence within the framework of PDF equa- 
tions focussing on deviations from Gaussianity. A joint analytical and numerical analysis shows 
that these deviations may be quantified studying correlations of dynamical quantities like pressure 
gradient, external forcing and energy dissipation with the velocity. A stationary solution for the 
PDF equation in terms of these quantities is presented, and the theory is validated with the help 
of direct numerical simulations indicating sub-Gaussian tails of the PDF. 



, Understanding the statistical properties of fully devel- 

L__Dped turbulence is one of the main challenges in classical 
rion-cquilibrium statistical mechanics [Il[2]. A major ob- 
stacle for the formulation of a successful theory starting 
fpom the basic equations of fluid mechanics is the gener- 
— . ically non-normal statistics of physical observables like, 

t" — velocity increments across scales. This non- normal 
statistic renders the applications of well-known theoret- 

fSj ical tools like renormalized perturbation expansions and 
the application of the renormalization group theory, which 
have been highly successful in the field of equilibrium 

' . I phase transitions and critical phenomena, rather difficult 
K*' [2 ■ Hence the development of a statistical theory based 

k>( on first principles remains one of the most fundamental 
problems of the field. 

_C^_' In contrast to the complex intermittent behavior of 
the velocity increment statistics, for which various phc- 
nomenological theories have been proposed (see, e.g., [1] 
for an overview), it is a well-known fact that the single- 
point velocity PDF in homogeneous isotropic turbulence 
is close to Gaussian [5], a fact which is often regarded 
as rather trivial. While Gaussianity often is assumed on a 
phenomenological basis, careful numerical and experimen- 
tal investigations, however, suggest slight deviations from 
Gaussianity |B]-[5] with a tendency to sub-Gaussian tails. 
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The question to be raised towards fundamental turbulence 
research consequently is, whether the PDF indeed deviates 
from normality and if so, why. A number of works aim 
at answering this question with a somewhat inconclusive 
outcome. Jimenez suggests sub-Gaussian tails taking into 
account the precise shape of the energy spectrum , while 
the instanton formalism is used by Falkovich and Lebedev 
|10| to argue in favor of sub-Gaussian tails depending on 
the external forcing. In contrast to that Gaussian PDFs 
have been obtained in the case of decaying turbulence by 
Ulinich and Lyubimov [TT] and later by Hosokawa [H] or 
in the case of the cross-independence hypothesis by Tat- 
sumi and coworkers |13j . The latter works make use of 
the statistical framework of the Lundgren-Monin-Novikov 
hierarchy jl4H16j with additional closure approximations. 

The scope of the present work is to give a comprehen- 
sive and clear answer to the above formulated question. 
Moreover, we shall demonstrate in the present Letter that 
already the nearly Gaussian behavior of the most funda- 
mental PDF is a consequence of the interplay of several 
non-trivial correlations. To this end, we make use of the 
LMN hierarchy combined with conditional averaging de- 
riving a PDF equation for the single-point velocity PDF. 
This allows to recast the closure problem of turbulence in 
terms of statistical correlations between the velocity and 
the various dynamical quantities like pressure gradient, 
external forcing and the rate of energy dissipation. While 
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this type of equation has also a long tradition in the con- 
text of PDF modeling of reacting flows (see, e.g., [ITfEO] ). 
a profound combined theoretical and numerical analysis is 
still lacking. This brings us to the structure of the present 
article. After introducing the PDF equation for homoge- 
neous isotropic turbulence, we will take into account sta- 
tistical symmetries to simplify its functional form. This 
will allow to derive an exact stationary solution, which 
expresses the velocity PDF in terms of the above correla- 
tions. It turns out that the unknown correlations, which 
show up in terms of conditional averages, are subject to a 
number of constraints, which motivate a Gaussian solution 
as a lowest-order approximation closely related to the as- 
sumption of statistical independence of the various quanti- 
ties. These theoretical aspects then will be supplemented 
by results from direct numerical simulations (DNS). The 
joint investigation then will show that deviations from this 
lowest-order approximation occur and that a subtle inter- 
play of correlations leads to a sub-Gaussian shape of the 
PDF. 

We start from the incompressible Navier-Stokes equa- 
tion 

d 

—u + u-\7u:^-yp + iyAu + F. (1) 
ot 

Here u{x,t) denotes the velocity field, p{x,t) is the pres- 
sure, 1^ denotes the kinematic viscosity and F{x,t) de- 
notes a large-scale forcing applied to the fluid to pro- 
duce a statistically stationary flow. PDF methods al- 
low to derive evolution equations for turbulent quantities 
right from the basic equation of motion. It can be shown 
[T4Hl6l[20l[2T| that the evolution equation for the single- 
point PDF f{v;t) — (^5{u{x,t) — v)) of the velocity for 
homogeneous turbulence may be derived taking the form 
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where the right-hand side is governed by the conditional 
averages of the pressure gradient, the external forcing and 
the so-called conditional dissipation tensor with respect to 
the velocity vector. Here we have introduced the sample 
space variable v for the velocity field. In the following we 
will restrict ourselves to the special case of homogeneous 
isotropic turbulence. In this case the PDF of the velocity 
vector is fully determined by the PDF of its magnitude, 
f{v), by the simple relation 



f{v) = ATTv^fiv), 



and eq. ([3]) may be recast taking the simple form 
H-|-$ 



dt dv 



(3) 



(4) 



Here the pressure and forcing term are expressed as 

Vp|t;) = H(u)'U n(u) tt-Vplw) 
{f\v) = <^{v)v $(i') =(m-F|u), 



i.e., as the conditionally averaged projections on the direc- 
tion of the velocity ii. The conditional dissipation tensor 
is characterized by two scalar functions 
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{v{k^uf\v) 
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which are related to the conditional average of the rate of 
energy dissipation, the enstrophy and the additional term 
(u{PJ' uf'\v^ related to the stretching and turning of the 
velocity vector by the velocity gradients Aij = . Due to 
the statistical symmetries eq. ^ represents an effectively 
one-dimensional problem with the appearing conditional 
averages depending only on the magnitude v of the veloc- 
ity vector, cf. eqs. (H])-®. The most interesting feature 
of eq. (m now is that it possesses a unique stationary 
solution 
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expressing the PDF of the magnitude of velocity in terms 
of the introduced conditional averages. Due to eq. ^ 
this already quantifies the PDF of the full vector. The 
statement of this exact result is that the explicit shape 
of the PDF depends on the correlations of the different 
dynamical quantities with the velocity on the single-point 
level. Of course, this result still contains unclosed terms; 
its value lies in the fact that the quantities determining the 
shape of the stationary PDF are identified and can now 
be discussed on the basis of physical arguments. Alterna- 
tively, the conditional averages can be expressed in terms 
of the two-point velocity PDF, such that the problem can 
also be formulated in terms of multi-point PDFs. This, 
for example, has been done in [T21II31I21] ■ We here rather 
follow a discussion in terms of the conditional averages. 
The precise knowledge of these functions now would allow 
to determine the exact shape of the velocity PDF. How- 
ever, they cannot be calculated explicitly without further 
closure assumptions, but numerical simulations or exper- 
iments may be used to obtain further insights. Neverthe- 
less, a number of constraints can be formulated for these 
unknown functions as upon integration the conditional av- 
erages have to reduce to ordinary averages. 



Av V H(u)/(w) 
Avv<^{v)f{v) 



{u • Vp) = 



(8) 



Av[\{v) + 2^i{v)]f{v) = iv 



u-F) = (e) (9) 
dui dui 







dxk dxk 
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It further can be shown on the basis of isotropy that 

dA I d£t I 

dv \v—i^ dv \v—{) 



(5) H(0) = $(0) = 0, A(0) = Ai(0) and 
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0, which further constrains the possible functional forms 
of the unclosed terms. 

The closure problem now prohibits an exact analytical 
treatment without further assumptions. Before we esti- 
mate the unclosed terms from our DNS data, we would like 
to demonstrate that Gaussian statistics can be obtained 
as the result of a number of simple modeling assumptions. 
The purpose of this short presentation, however, is not 
to claim the velocity statistics to be exactly Gaussian. 
It rather shows that physically intuitive, but in the end 
inappropriate, assumptions can lead to quite reasonable 
results. 

To proceed with an analytical treatment of the PDF 
equation, one now, on the one hand, could assume de- 
coupling of the velocity field that varies on large scales 
and the dissipation field, which varies on much shorter 
scales. As a consequence, the conditional dissipation ten- 
sor is a constant isotropic tensor. The pressure gradient 
cannot be treated with such a simple argument, because 
it is known to be a strongly fluctuating, but long-range 
correlated quantity. Still, to comply with the integral con- 
straint ([5]) , the most simple modeling assumption is a van- 
ishing contribution, which is related to the assumption of 
statistical independence of the pressure gradient and ve- 
locity field. Finally, the external forcing and the velocity 
field both vary on large scales, such that mutual depen- 
dences have to be expected. These considerations lead to 
a lowest-order approximation (in terms of powers of v) 
compliant with the above constraints, taking the form 

Iloiv) = (11) 

Mv) - (12) 

Xo{v)^M = ^, (13) 

where cr = Umis = \f^^ denotes the standard deviation 
of the velocity. That means, in this lowest-order approxi- 
mation the conditional pressure gradient vanishes, and the 
external forcing depends linearly on the velocity on the 
statistical average. The terms related to the conditional 
dissipation tensor are constant, viz. statistically indepen- 
dent of the velocity field. Insertion into eq. ([7]) yields the 
Gaussian solution 

1. e., under a lowest-order ansatz including the assumption 
of negligible correlations we have found a closure of cq. 

that predicts Gaussian velocity distributions0 In turn, 
the existence of non-negligible correlations will generically 

-"^In fact, any closure with Ilo{v) + ^q{v) = and Ao{'u) = 

MoC") = leads to a Gaussian velocity PDF. The present one is 
additionally consistent with the integral constraints lISll- lllOII and the 
corresponding energy budget equation for stationary homogeneous 
isotropic turbulence. 





i?A 


L 




512^ 


112 


1.55 


2.86 0.543 10-^ 0.103 2.03 



Table 1: Major simulation parameters. Number of grid points 
A'''', Reynolds number based on the Taylor micro-scale Rx, in- 
tegral length scale L, large-eddy turnover time T, root-mean- 
square velocity itrms, kinematic viscosity i/, average rate of en- 
ergy dissipation (e). femax?? characterizes the spatial resolution 
of the smallest scales. 

cause deviations from Gaussianity, which will be discussed 
in the following. In fact, the conditional averages obtained 
from our numerical simulations will significantly deviate 
from the modeled ones above and together will lead to 
deviations from Gaussianity. 

The numerical simulations are obtained with a standard 
pseudospectral code in a cubic box with periodic boundary 
conditions, for which the main simulation parameters are 
summed up in table [1] The results were checked within a 
range of Reynolds numbers up to R\ « 220 at a resolution 
of 1024'^ and for several large-scale forcing methods. The 
results of the analysis proved to be robust in the sense that 
all of the presented results do not vary strongly within the 
range of investigated Reynolds number and the different 
types of forcing schemes. It would, however, be of gen- 
eral interest to obtain the statistical quantities presented 
in the following for fiows at significantly higher Reynolds 
number in order to study Reynolds number effects. Due 
to its excellent statistical quality we present the results 
of a 512'^ simulation, where the statistics has been evalu- 
ated for a run longer than 150 large-eddy turnover times 
sampled sufficiently often (about 200 snapshots have been 
stored) to form a proper statistical ensemble. As deriva- 
tives of the velocity field like the energy dissipation have 
to be evaluated, special emphasis has been put on an ap- 
propriate resolution of the smallest scales of the flow. 

In the following we will evaluate the terms that de- 
termine the shape of the velocity PDF according to eq. 
([7]). The conditional pressure gradient term and the con- 
ditional forcing term are shown in flg. [T] The most striking 
feature is that the pressure contribution is non- vanishing. 
For low values it is negative indicating that fluid parti- 
cles on average will be decelerated due to the (negative) 
pressure gradient. For high values of velocity, however, the 
term becomes strongly positive, which correspondingly in- 
dicates an acceleration of fluid particles with high veloc- 
ities due to this term. If the pressure gradient term is 
non-vanishing, the observed zero-crossing is necessary to 
fulfill the integral constraint ([5]). These findings show that 
pressure contributions may not be neglected as opposed 
to the analysis presented in [12l[T3]. The forcing term is 
positively correlated with the velocity and hence also has 
an accelerating effect on a fluid particle on the statistical 
average. For low values of velocity this effect is approxi- 
mately linear and seems to saturate for higher velocities, 
which suggests that the forcing has a minor effect on the 
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Figure 1: The conditional averages n(i;), A(i;), and the 

PDF f{v) shown for reference. The conditional averages show 
an explicit u-dependence indicating statistical correlations. 



Figure 2: The conditional averages ^(i;), \{v) and the PDF 
f{v). The statistical correlations increase in the tail of the 
PDF. 



tails of the PDF. It can be shown in the framework of 
PDF equations that the conditional acceleration, i.e. the 
conditional right-hand side of the Navier-Stokes equation 
([T]), vanishes [23]. For isotropic turbulence this implies 

n(w) + + A(i;) = 0, (15) 

by which wc have introduced the viscous forces A(i') = 
{vu ■ Au\v). This term indeed has a decelerating effect 
on a fluid clement as can be seen from fig. [T] On the 
whole the pressure gradient and the external forcing are 
balanced by viscous forces. 

For the stationary solution of the PDF equation we ad- 
ditionally have to take into account the terms from the 
conditional dissipation tensor. For PDF modeling and in 
the case of some closure theories [TIHT^ITTlfT^ this ten- 
sor is taken to be a constant isotropic tensor of rank two, 
compliant with our lowest-order Gaussian approximation. 
Fig. [2]shows the outcome of our DNS analysis. Indeed, we 
find the approximate relation A « /i for the two eigenval- 
ues holding especially well for moderate values of velocity 
indicating that this tensor is approximately isotropic. It is 
furthermore nearly constant for low values of velocity, such 
that the above approximation is appropriate for this region 
of the PDF. For high velocities, however, stronger statisti- 
cal correlations are apparent showing that high velocities 
on average are correlated with dissipative events. Conse- 
quently, the tails of the PDF will be influenced by these 
correlations leading to a more rapid decay (recall, e.g., 
that A appears in eq. ([7]) in the denominator). The fact 
that the conditional dissipation tensor is approximately 
isotropic indicates the absence of directional correlations 
of, e.g., the energy dissipation field with the velocity field. 
In view of the aforementioned results the assumption of 
a decoupling of large- and small-scale fields may be re- 
garded as a rough approximation leading to a Gaussian 
closure, and one is likely to neglect important features of 
the statistics of turbulence under this assumption. 



These results can be used to motivate improved closures 
by modeling the unclosed terms, such that they resemble 
the functional shape observed in the DNS. This, of course, 
brings up the question if the observed results are of uni- 
versal nature, and it would be interesting to perform the 
analysis also for different flows apart from the idealized sit- 
uation of homogeneous isotropic turbulence. In the sense 
that we expect statistical correlations between, e.g., the 
dissipation and velocity field also in more complex situ- 
ations, we think that the obtained insights are of more 
general nature. 

Now having obtained all the unclosed terms from the 
DNS, we can explicitly evaluate eq. Q, which is pre- 
sented in fig. [3l As all of the presented results leading to 
eq. ([7]) are exact, the PDF calculated from eq. ([7]) per- 
fectly collapses with the velocity PDF directly determined 
from the DNS data. This consistency check proves that 
the interplay of statistical correlations of the velocity with 
the dynamical effects of the pressure gradient, external 
forcing and dissipative terms causes the deviations from 
Gaussianity and especially the sub-Gaussian tails of the 
PDF. In view of the pronounced statistical correlations 
for large values of velocity it is actually surprising that 
deviations from Gaussianity are comparably small, which 
shall be explained in the following. A closer look at eq. 
([7]) shows that deviations from Gaussianity may arise for 
different reasons. Considering the exponential factor, de- 
viations from Gaussianity appear if ^ is non-constant and 

is non-linear. These two quantities are presented in 
fig. U along with the corresponding quantities from the 
analytical Gaussian closure approximation (fTT|) - P^ . It 
turns out that the strong correlations observed in figs. [T] 
and [2] tend to cancel to a certain extent, such that the 
quotient of the sum of the (negative) pressure gradient 
and the external forcing over the eigenvalue A of the con- 
ditional dissipation tensor is approximately linear and j 
is approximately constant. The former is possible since 
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Figure 3: Comparison of the PDF /(n) with the evaluation of 
eq. (O- The reconstructed PDF collapses with the directly 
estimated one. The excellent agreement demonstrates the con- 
sistency of the theoretical results. A comparison with an angle 
integrated Gaussian shows that, although being close to Gaus- 
sian, significant deviations exist. 

regions of high velocities are on average associated both 
with large accelerating forces due to the pressure gradient 
and strong dissipative events; in these regions the forcing 
plays a minor role compared to the other dynamical influ- 
ences. The latter approximate behavior comes due the fact 
that only weak directional correlations between the veloc- 
ity field and the dissipation field occur. The slight devi- 
ations from linearity and constancy, respectively, of these 
quantities with their particular functional form eventu- 
ally contribute to deviations from Gaussianity. Additional 
deviations arise due to the functional shape of A, which 
appears in the prcfactor of cq. ([7]), as the Gaussian clo- 
sure approximation p3p here suggests a constant function. 
By and large the subtle interplay of the pronounced ob- 
served statistical correlations causes rather moderate but 
still significant deviations from Gaussianity towards sub- 
Gaussian tails. 

To conclude, we have analyzed the single-point statistics 
of fully developed turbulence in the framework of the LMN 
hierarchy. Taking into account statistical symmetries, we 
have derived an expression for the stationary PDF, which 
highlights the connection between the shape of the PDF 
and correlations of the velocity with dynamical effects. 
This stationary solution is an exact result derived directly 
from the Navicr-Stokes equation. One of the strengths of 
this approach is that the closure problem of turbulence 
is expressed in terms of correlations on the single-point 
level rather than the traditional approach in terms of a 
coupling to higher-order moments or multi-point PDFs, 
allowing for a physical discussion of the unclosed terms. 
Furthermore, we have found Gaussian solutions applying 
simple physical closure assumptions consistent with the 
constraints on the conditional averages. 

To explain the observed deviations from Gaussianity, we 
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Figure 4: The quotients (H + $)/A and ^/A appearing in eq. 
(O show an approximately linear and constant functional form, 
respectively. The corresponding quotients expected from the 
analytical Gaussian closure approximation ([TT])-([T3| are shown 
for comparison. 

have used DNS data to study the conditional averages that 
contribute to the shape of the PDF finding non-negligible 
correlations. Evaluation of the stationary solution of the 
PDF equation proved to be in perfect agreement with the 
directly estimated DNS data. The results demonstrate 
how a combined theoretical and numerical investigation 
can be used to give insights into the details of the statistics 
of turbulent flows. Apart from giving a comprehensive pic- 
ture of the single-point velocity statistics, the results may 
motivate improved closure assumptions for PDF modeling 
by taking into account the functional shape of the condi- 
tional averages observed in DNS. It remains an interesting 
question for further research if the observed correlations 
carry over to more complex flow situations like, e.g., wall- 
bounded turbulence. 

The theory can also be applied to decaying turbulence. 
An extended analysis including more details on the theo- 
retical derivation, application to decaying turbulence and 
a comparison to the vorticity statistics will be the topic of 
a detailed report published elsewhere (24) . 

^ ^ ^ 
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